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1 Introduction 



This paper is devoted to the study of Bethe Ansatz equations (BAE) for the Isotropie spin-2 
Heisenberg chain. The model describes interacting spins, situated on the sites of a periodic lattice 
A'^. There has been a profound development in physical interest on the thermodynamical properties 
of their Solutions. In recent years a strong increase of mathematical attention on this subject arises 
from the study of integrable Systems via quantum inverse scattering method. In most cases the 
results obtained are mainly concerned with the lattice size N tending to infinity. However the 
finite-size problems, or the finite-size corrections, should be of mathematical interest, as well as 
essential to the understanding of BAE, which is the basis for all exact Solutions of electronic models 
in one dimension. We intend to investigate here the finite-size BAE for Heisenberg XXX model in 
a rigorous mathematical way, which will serve as a Standard theory for our subsequent works on 
quantum integrable Systems. 

In his remarkable works Baxter discovered a link between the quantum spin models and 
transfer matrices of 2-dimensional Statistical lattice models. The quantum inverse scattering 
method ^ provides a natural understanding of the role of Bathe Ansatz in the problem of spec- 
trum of model Hamiltonians. For Heisenberg XXX Hamiltonian, Bethe vectors have some special 
characteristic properties: they are eigenvectors of of the corresponding transfer matrices with poly- 
nomials as the eigenvalues, and also the highest weight vectors for the global s/2-symmetry of the 
Hamiltonian. In the investigation of Solutions of BAE in the large N limit, there is one fundamen- 
tal assumption, the so-called string hypothesis (Takahashi Q). We observe that a certain feature 
of these string structures can be understood through BAE, while the counting of states Q based 
on this string hypothesis remains still valid for the 2-particle sector of any finite System N. The 
analytical Solutions for sectors other than two particles are hard to obtain except a very small size, 
e.g. = 6. However for the ground state, one can study the problem by analysing the corre- 
sponding logarithmic form equation, where the fixed point theory can be used for the existence of 
Solutions. The uniqueness of the ground state Solution should be mathematically expected, but 
only strong indication can be obtained here. Such program is now under progress and partial results 
are promising. 

The Organization of this paper is as follows. In Section 2, we recall necessary definitions in 
the theory of quantum integrable Systems Q and give various characterizations of Bethe roots 
for Heisenberg XXX model. In Section 3, Bethe roots as A^ — > 00 are discussed, and the string 
structure for a large A^ is derived from BAE. In Section 4, we study some problems on Heisenberg 
XXX model of a finite lattice size. Here the BAE for the two-particle sector is solved analytically, 
and mathematical structures of BAE on the ground state in antiferromagnetic case, as well as the 
equivalent equation for its logarithmic form, are discussed. We establish rigorously the existence 
of ground state via the fixed point theory, and also the uniqueness of the Solution for a small N. 
In Section 5, we present an illuminating ( but not a mathematically rigorous ) argument on the 
uniqueness of ground state for a large finite System. We have written this note in a mathematical 
style, and hope that in process it would not be much difficult to read for both mathematicians and 
theoretical physicians. 

2 Characterizations of Bethe States 

The Hamiltonian of Heisenberg XXX model is given by 

•^^11 22 33 

-f^XXX = + «"nC^n+l + (^n(^l+l " 1) > J € IR . 

n=l 
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Here er-' 's are Pauli matrices. In this paper, we shall always assume the size to be even with 
the periodic boundary condition ( + 1 = 1) . The Operator Hxxx acts on the Hilbert space of 
physical states Hat, 

N 

= <^ hn , K := (C^ for all n. 

n=l 

The link between the above quantum XXX System and a 2-dimensional Statistical lattice model is 
described as follows ||5|. Dehne the Operators of (D^ /i 



with the first factor (D^ as the auxiliary space, and the second factor h as the quantum space with 
the basis 











, i->H 





1+ >= 

The R{X) satisfy the following Yang-Baxter equation: 

7^(A - ij){R{X) R{fi)) = {R{fi) R{X))TZ{X - jj), 
where 7^(A) is the 4x4 numerical matrix defined by 



TZiX) 



f 1 \ 
b(A) c(A) 
c(A) b(A) 



, b(A) = ^, c(A) = ^ 

\ 1 J 



Using R{X), we introduce the local transform matrix 

L„(A) = 1 • • • (2) R{X)nth • • • O 1 
as the Operators of (C^ Hjv having R{X) at the n-th side. Define the monodromy matrix 

Fm{x) = Lm{x)Lm.,{x)---l,{x) = ^j^j ^1^1 ^ . 

Then we have 

7^(A - ^i){FN{X) ® FNifi)) = (FNiß) F7v(A))7^(A - fi). 

The matrix entry elements A{X), B(X),C{X), D(X) of the monodromy -F(A) are Operators of Hat, 
which generate the so called ABCD algebra: 

AiX)A{^i) = A{i^)AiX), B{X)Bif,) = B{f,)B{X), 

CiX)Cif,) = C(m)C(A), DiX)Dip) = Dif,)DiX), 

A{ß)B{X) = f^,xB{X)A{fi) + gy^B{fi)A{X), A^B, 

D{X)B{fi) = f^,xB{fi)D{X) + gx,^,B{X)D{fi), B ^ D, 

C(A)A(^) = /^,a^(m)C(A) + 5A,^^(A)C(/i), A^C ^> 

Cif,)D{X) = f\,,xD{X)Cif,) + gx,,Ml^)C{X), C^D 

C{X)B{j,) - Bif,)C{X) = gx,,,iA{X)D{j,) - A{i,)D{X)) = g,^x{D{X)A{^,) - D{i,)AiX)), 

D{X)A{ti) - A{t,)D{X) = gxAB{X)C{f,) - i?(^)C(A)) = g^,x{C{X)B{t,) - C(/.)i?(A)). 
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where 

1 b(A-^) 
^"'^ c(A-m)' ^^''^ c(A-/x)' 
Taking the trace of the monodromy, one obtains the transfer matrices 

r^(A) := trFjv(A) = ^(A) + D{X) , 

which form a family of commuting Operators of Hjv: 

rAr(A)-r^(^) = rAr(//)-rjv(A). 

The Hamiltonian of XXX model is related to the transfer matrices by 

//xxx = -^(iJ^logr(A)U^. -iV) . 



Define the pseudo-vacuum 



N 

'■= Vn G Hjv , Vn = \+ > for all n. 

n=l 



Then we have 



A{\)nN = ix + ^)^nN , L>(A)Ojv = (A - i)^j^iv , c{x)nN = o, 

For complex numbers < j < l, we consider the vector 



*iv(Ai,...,AO := l[B{Xm)nN e Hat, 



A(A; Ai, . . . , Xi) = A(A; Ai, . . . , A^) , A(-A; -Ai, . . . , -A,) = A(A; Ai, . . . , A,) , 

hence 



{A,}, CR ^ A(A;A,) = A(A;A,) 

{A,}i = {-A,},- ^ A(-A;A,) = A(A;A,-: 



By the relations 



nA Xm ^ _ 1 1 \ ^ '^m TT A Aj,^ ^ _ -I I \ ^ ^5 



A-A^ ^=iA-A^ A-A^ ^=iA-A^ 



where 

dm = n ~ \ — ZT") ' = ^ n + 



(2) 



m=l 

and define the function of A, 

A(A; A„ . . . , A,) := (A + \f ^ + (A - II ^\\^' " (3) 

m=l m=l 

One has 
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one has 

A(A; Ai, . . . , = (A + 1)^ + (A - '-f + ^ -rniX + + b^{X - 

L L Ä — Am 

m=l 

The criterion for A(A; Ai, . . . , A;) to be an entire function of A is now equivalent to A^'s satisfying 
the Bethe Ansatz equation ( BAE ) : 

(^)^= n ^^^^ , X,e(E, j = l,...,l. (5) 

~ 2 m=l,m^j ^1 * 

In this Situation, the difference of A(A; Ai, . . . , A^) and (A + + (A — is a polynomial of degree 
at most — 1. Note that (A + + (A — is the eigenvalue of the transfer matrices Tj\[{X) for 
the pseudo-vacuum , 

T^(A)17,v = ((A + + (A - '-f)QN ■ 

The variable A can be interpreted as the "rapidity" of a "particle" with its momentum p{X) defined 
by 

jp{\) _ ^ + f 

or equivalently , 

The scattering amplitude S{X, fi) of two particles is 

fi — X — i 

with the scattering angle satisfying the relation 

2cot^^^ = cot^-cot^. 
2 2 2 

Now BAE (^) also takes the form 

i 

^ip{\k)N = — Y[ e*'^^''**'^'"-* , k = 1, . . . ,1. 

m=l 

The right hand side is the /-particle scattering amplitude in terms of two-particle ones, a fact that 
manifests the integrability of the model. By the relation 

S(A,/.)5(/.,A) = l, 

one obtains 

k=l k=l - 

The vector ^'^^(Ai, . . . , A;) is called a Bethe vector when Aj's form a Solution of BAE. 
Denote 

1 ^ 

Sj = isT.< (j = l,2,3); S± = Si±iS2. 

^ n=l 
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S±,S3 form a basis of 5/2 ((C) acting on Htv- One has 

= i E/=l[ö-^ O-^ ® crl]aux = ^[(^^ , Ln{X)]aux 

here the Lie-product [M, 6] of a 2 x 2 matrix M, 

\ M21 M22 ) 

with an element 6 in an arbitrary ring is defined to be the 2x2 matrix 



[MM 



[Mn,6] [Mi2,b] 



^ [M21M W22A ) ' 
and the Lie-bracket of the last term is on the auxiliary space. It foUows 

[FN{X),Sj] = ^[a^,FN{X)]aux 



and 



[A{X),S3] = [Di\),Ss] = [C{X),S+] = [C(A),5_] = 0, 
[BiX),S3] = BiX) , [CiX), Ss] = -C(A), 
[B{X),S+] = -[C(A),Ä_] = -A{X)+DiX), 
[A{X),S+]=-[D{X),S+] = C{X) 
[A{X),S^] = -[D{X),S^] = -B{X) 



Therefore 



[TN{X),Sj]=0 , 

which States the sZ2-symmetry property of the Hamiltonian i?xxx- It is easy to see that 

N 

S+Ün = , S3ÜN = -^^N 
and the spin of "^NiXi, • • • , A;) is ^ ~ ^ > 

Ss^NiXi, . . • , AO = (y - /)*iv(Ai, • • • , Xi). 

The foUowing conditions for a Bethe vector should be well-known speciahsts, but we could not find 
exphcit references. Here we give the details of the proof. 

Proposition 1 . ^]v(Ai, . . . , A;) is a Bethe vector if and only if one of the foUowing equivalent 

conditions holds : 

(i) {Aj}^-=i satisfies BAE. 

(ii) The function A(A; Ai, . . . , A;) has no pole at a finite value of A. 

(iii) ^'Ar(Ai, . . . ,Xi) is a common cigcnvcctor of transfer matrices Ti^{X). 

(iv) ^'7v(Ai, . . . , A;) is a highest wcight vector ( of spin ^ " ^ ) ^01 the si2-i'epresentation. 
Proof. The equivalence of (i) and (ii) is known before. 
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(i) (iii). By the relations @) one can move A{X) and D{X) through B{Xm) to Oat in the 
expression 

i l 
Tn{X) n Bi^m)^N = {A{X)+DiX)) n BiXm)nN . 

m=l m=l 

The resulting form becomes a combination of 

i 

A(A;Ai,...,A;) n BiXm)nN 

m=l 

with the forms 

l 

Afc(A;Ai,...,A05(A) [] ^(^'n)^^w , (A: = 1,...,0. 

The expression of A(A; Ai, . . . , A;) is given by dH) by taking account only the first terms on the 
right hand side of the commutation relations of A,B and D,B in (|l|). If we use second terms of 
(H) on the commutation relations of A{X) , D{X) with B{Xi) , and then first terms of (||) on the 
commutation relations of ^(Ai) , -D(Ai) with B{Xm) for m > 2, by the relation gx^^ = —g^^x , we 
obtain the expression of Ai(A; Ai, . . . , A;): 

Ai(A; Ai, . . . , A,) = 5a„a[(Ai + "-f {[ /a,a„ - (Ai - '-f {[ /a^.aJ • 

m=2 m=2 

Since i?(Am)'s are commuting Operators, by a suitable permutation of the indices m, one concludes 
the expression of Afc(A; Ai, . . . , A;) which is given by 

i ■ i 

Afe(A;Ai,...,A0 = <7A„A[(A,. + ^)^ U fx,,x^-{Xk-\f J] /a^aJ • 

Therefore ^'Ar(Ai, . . . , A/) is an eigenvector of T]\f{X) provided {XkYi^^i satisfies the relations 

Afc(A; Ai, . . . , A;) = , for k = l,...,l, and all A, 

which is equivalent to BAE. Therefore we obtain the result. 

(i) <-> (iv). Using the relations between B{X) and Sj, one has 

i 

S+^n{Xi, ■■■,Xi) = T. B{M) ■ ■ ■ ß(A,-i)(A(A,) - D{Xj))B{Xj+i) ■ ■ ■ B{Xi)nN. 

By using (|l|), S'+^'iv(Ai, • • • , A;) can be expressed as a combination of forms 

Mfc(Ai, . . . , Xi)B{Xi) ■ ■ ■ i?(Afc_i)ß(Afc+i) • • • B{Xi)nN , {k = l, . . . ,1). 

Taking account only of first terms on the right hand side of the commutation relations between 
A,B and D,B in dTj), one obtains 

Mi(Ai, . . . , Az) = (Ai + -)^/Ai,A2/Ai,A3 • • • /Ai,A, - (Al - ■^)^ fx2,XifX:i,Xi ■ ■ ■ fXi,Xi- 

Since i?(Am)'s are commuting Operators, by the symmetry argument, one also has 

i ■ i 

Mk{Xu...,xi) = iXk+'-f n /a.,a™ - (Afc - n /a™,a,. 
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The vanishing of Mk{Xi, . . . , A;) is now equivalent to BAE. Hence we obtain the equivalence between 
(i) and (iv). □ 

Remark. The Bethe states are in an one-one correspondence with the Solutions of BAE. Indeed 
one has the following equivalent statements for Bethe vectors: 

'I'Ar(Ai, • • • , Az) = s-^NiK^ • • • ; K) some s G (D - {0} 
^ A(A;Ai,...,AO=A(A;Al,...,A',) 
<^ l = k and {Ai, . . . , A^} = {A'i, . . . , A'^,} . 

It is obvious that two linear dependent Bethe vectors ^'Ar(Ai, • • • , A^) and ^'Ar(A'i, • • • , A'^) have the 
same eigenvalues: 

A(A;Ai,...,AO=A(A;A;,...,A'fc), 

and by (13), this means 



a^(A + §)^ + b^(A-§)^ ^ ^ a;.(A + §)^ + 5;.(A- 
A — A777 . A — Aj 



2' 



m=l j=l J 



where 



„I 



Note that by BAE, 



Then we have 



&m = <;=^ Am = -ö^ =^ O-m 7^ . 



(A + ^nET^-ET^) = -(A-V(ET^-E J 



" -1 A — Am 1 A — AI- 2 A — Am ^ A — AI- 

m=l j=l J m=l j=l J 

By multipling nL=i(A — Am) 11^=1 (A — A^-) on the above equation and from l + k < N, we have 

l k l l 1 k u 

_ ^ "j _ Q ^ _ ^ j _ Q 

-1 A — Am -^A — A,- iA — Am ^ A — A„- 

m=l j=l J m=l j=l J 

This implies that each Am is equal to some X'j. Similarly A^ is equal to Am for some m, hence 

l = k and {Ai, . . . , A;} = {A'i, . . . , A'^} . 
So we obtain the conclusions. □ 



3 String Hypothesis of Roots of Bethe Ansatz Equation 

In the description of Solutions of BAE as the size tends to infinity, one assumes the string 
hypothesis which claims any Solution consists of a series of strings in the form 

A = X + - fc) + 0(e-'3^) k = l...n xGM, 

for some ß > 0. In this section we shall discuss the relation between the string structure and BAE. 
First we dehne certain notions needed for our purpose. A complex number z is called an asymptotic 
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limit point of a sequence of Bethe roots {A^^^}'-^^-' if for some choice of j{N), the following relation 



holds: 



A^.JJ^ = z + 0(e-^^) as ^ oo . 
For the convenience, we shall simply write 

Xj^^z as N —>■ oo 

when no confusion could arise. When the asymptotic hmit point z is non-rcal, i.e. Im(2:) 7^ 0, the 
above element A^.^^ is called a complex root in the Bethe Solution {Xj^^'j^j^^ . Moreover, if there 

exists another sequence A^^^ whose asymptotic hmit point is the real part of a complex root, A^^^ 
is also defined to be a complex root. 

A more precise description for string hypothesis states as foUows: A Bethe Solution for a large 
N always lies in a collection of Bethe Solutions {Xj^^Yj^^, ( iV 3> 0), such that every root belongs 
to a suitable convergent sequence of the form: 



X) ' ^ X + i{ — k) as N ^ 00 . 



3 

The collection 

Tl — 1 

x + i{^^-j) , 0<j<n-l xGR, (7) 

will be called a string of length n with center x. We shall derive this string structure of Bethe roots 
as N tends to infinty under the following additional ( somewhat unpleasant ) condition: 

Hypothsis (H) : " There exists a positive integer Nq such that any Bethe root {Aj}j-^i for 

a large size always lies in a sequence of Bethe Solutions {X^^^Yj^^ ,N > Nq, which tends to an 

asymptotic configuation consisting of Clements in real axis with a finite number of complex numbers. 
Futhermore the number of complex roots of Bethe Solutions remains constant in the limit process 
as AT ^ 00." 

First, let us consider the simplest case for Z = 1. BAE is given by 

(^)-^l.Ae(D. (8) 

As N tends infinity, the above equation becomes 

which simply means the momentum p{X) being real, i.e. the real rapidity A, and the Bethe State 
^n{X) is called a magnon state with the energy given by ^^^/^i^ . 

For / > 0, we describe the following lemma which somewhat suggests the conjujate Symmetrie 
nature of a string. 

Lemma 1. Let {XjY^^i be a Solution of BAE with iV — 00 . Assume 



x + i{y- j) j <n 
GM j > n 



for some positive integer n and real numbers x, y . Then Ai, . . . , A„ form a string of length n with 
center at x. 
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Proof. We have 

i ^ I AI.. 1 



f{h±l= -+;(y-^)^ 1^1 = 1 for.>n. 

7=1 -^i - § X + i(y - n - i) Afc - § 



By and letting N ^ oo, one has 



2> 



x + i(y - n - i) 

which implies y = Therefore {Aj}"^]^ is a string of length n. □ 
For 1 = 2, BAE is given by 

A]_-M/2 yv _ Ai - A2 + i A2 + V2 jv _ A2 - Ai + i 
^Ai-V2 ~Ai-A2-i ' ^A2-i/2^ ~ A2 - Ai - i ' 

We have 

Ai + i/2 yv A2 + V2 jv ^ 
^Ai-i/2^ ^A2-i/2^ 

Either both Aj's are real, or both not. If Ai and A2 are real, \I'Ar(Ai, A2) is called a 2 magnon state 
and its energy equals to the sum of those of 1 magnon states ^'Ar(Ai) and \I'Ar(A2). In the case for 
both Ai, A2 not real, we may assume 1^^-377! I is greater than one, hence liniTv^oo I aI-^s ~ 
The first relation of BAE implies Ai = A2 + i- By Lemma 1 , Ai and A2 form a string of length 2, 
hence Ai = A2 = x + | for some x G IR. In this case ^'Ar(Ai, A2) is called a bounded state, and its 
energy is equal to ^j^^^^ • 

For l = 3, let {Ai, A2, A3} be a Solution of BAE. By the realtion 

j=i 2 

either all the Aj's are real, or at least two of them are not real numbers. The former case is the 3 
magnon state 'I'Ar(Ai, A2, A3). Otherwise we may assume Ai and A2 are not real with 

I Ai + ^ I ^ , I A2 + f 



Ai - § 

By BAE for Ai, A2 with N ^ 00, we have 



>1, |-^|<1. 



J (Ai-A2-i)(Ai-A3-i) =0, 
I (A2-Ai + z)(A2-A3 + i) =0 . 

Either Ai = A2 + i, or A3 = Ai — i = A2 + i- For A3 G H , by Lemma 1, either {Ai, A2} forms a 

A3+— 

string of length 2, or {A2, A3, Ai} forms a string of length three. For A3 IR, |t — f\ 7^ 1. By BAE 

A3-2 

for A3 with large limit, one has 

{Xi, X2, X3} = {z — i, z, z + i} for some z S (C, 

hence it is a string of length three by Lemma 1, which contradicts Xj IR for all j. In this way we 
have determined the structure of {Xj}j^i as tends infinity, which is composed of either 3 real 
roots, a real with a 2-string , or a string of length 3 . 



10 



For / > 4, the procedure we employ above is not sufficient to derive the string structure of Bethe 
Solutions, e.g. one needs to exclude a chain like 

^2 % 1 3z 

Ai = — , A2 = — , A3 = - , A4 = j , 

in the l = 4 Solutions. Now let {XjYj^i be a set of Bethe roots. Assume 

fj-k ■■= fJ- - {k - l)i £ {XjYj^-j^ , k = l,...,m (9) 
where ^ is a complex number and m is a positive integer greater than 1. Denote 

Mtl = {A,}^=i - MT=i ■ 

By multiplying the equations for fif^ in BAE, we have 

_i_ i m l—ni . ■ 

/ /^l "T 2 NiV = TT TT t^k-^h + i i^-^Q^ 



From 



one has 



., i (2m— l)i 



1/^1 + ii 1 1 1 T/\ m-l m-1 m - 1 

I tI = 1' >1' <1 ^ Im(^) = ^— , >^^, 

/^m ~ 2 Z Z Z 



respectively. Since the right hand side of (|10D is equal to 



from our assumption (H) one has 

Im(/i) > ^^^^^ =^ limAr_»oo T^h = A^m " ^ for some h 
Im(/x) < ^^^^^ liniAT^oo i^h = fJ-i + i for some /i , 

Im(/i) = ^^^Y^ =^ /ii, . . . , /im form a string of length m . 

As a consequence, if {/Ui, . . . ,/im} is a subcollection of {Ajl^-^^ which is maximal among chains of 
the form (|^, it must be a string of length m. Hence {XjYj^i is an union of reals roots and a finite 
number of strings of length greater than one in large limit. Therefore we obtain the following 
conclusion. 

Proposition 2 . Let {XjYj^i be a set of Bethe roots for site A^. As A^ — > 00, {XjYj^i is composed 
of real roots toghter with a finite number of strings with length greater than one. □ 
Remark. The centers of strings in a set of Bethe roots are indeed all distinct, i.e. Pauli principle 
holds for Bethe vectors, for the argument see e.g. 0] . 
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4 Bethe Ansatz Equation for a Finite Site A'^ 

In this section, we discuss the Bethe structure for a finite size A'^. 

Example 1. Betlie vectors for 1 = 1. BAE (|8|) is described the relation: 

^^^=ü;\ l<A;<iV-l, CO := e'^'/'' , 
X-i/2 ' - - 

hence 

A = cot— , 1< k < N - 1 . 

N ~ ~ 

By (H) , we have 

B{X)nN = i(A - '-)^-^ ^''"(1 ^■■■^^nth^---^ ^)^N , for A = cot ^ . 

71=1 

Hence the Bethe state ^Ar(cot ^) is a multiple of the vector 

N 

Uj''"'\+ >(g) ■ ■ ■ \+ >(g)\- >nth (S>\+ > • • • 1+ > 

n=l 

for 1 < k < N — 1, and all these Bethe 1-vectors form a subspace of Htv of dimension — 1. Note 
that the above vector for A; = corresponds to A = oo, which is the Solution of jr^yl = 1- n 
Example 2. Bethe Solutions for 1 = 2. Consider the change of variables: 

- = ^, A = l^. (11) 



We have 
and 

The / = 2 BAE becomes 
which is equivalent to 



A-| 2z-l 



1_A + | 



-z A-§ 



A e KU {oo} \z\ = l ■ 



N_ ZlZ2 - 2Z1 + 1 NN_. 
^ ZIZ2 -2Z2 + 1 ^ ^ 



f ziZ2 =uj^ , cu = e^^*/^ , < A; < iV - 1, 

] ^^-1 = i^'-'-2zi+l 

Note that if (21,22) is a Solution of the above /c-th equation, so are {z2,zi), and (=,=). Hence 
both zi, Z2, = , = are Solutions of the equation 

{üü'' + l)z^ - 2uj^z^-^ - 2z + (w*^ + 1) = 
for some k. Now we are going to determine its Solutions. Set 

e'P = üj-'^l^z , 
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then 



The above equation becomes 



COS ^e*^" - (-l)^e*(^-i)P - e'P + cos ^ = 0, 

For k odd or 0, p = 0, tt satisfy the above equation, and they correspond to the Solutions of BAE 
with zi = Z2 , hence not in our consideration. Note that the function fk{p) depends only on the 
parity of k, which can also be expressed by 



fkip) 



-'"»^-^;^(^^ for odd k 

sm{Np/2) 

cos((jy-2)p/2) £ even k 

cos (Afp/ 2) 



Hence fk{p) has the period 2tt with the symmetries fk{—p) = fk{p)-,fk{j^ ~ P) = ~fk{p)- So it 
suffices to consider the Solutions with < Re(/o) < vr and Im(yo) > . First let us determine the 
real Solutions of ( |T^ ) for < p < vr. Claim: /^(p) > 0. 
Since 

' ^-^-^^^[(iV_l)sinp-sin(iV-l)p] for odd k 
2cosHNp/2) [(^ - 1) P + ^^^(^ - 1)^] *5ven k , 



f'kip) 



it suffices to consider the region of p with 

1 

sm p < j^—j , < p <TT , 

which implies 

By fk{TT — p) = —fk{p), it needs only to consider the region 

hence /^(p) > for even k. For odd k, we have 

-^[(iV-l)sinp-sin(A^-l)p] = (A^ - l)[cosp - cos(A^ - l)p] >0 
dp 

hence 

(iV-l)sinp-sin(iV-l)p > 0, 

which implies /^(p) > 0. Therefore fk{p) is an increasing function on each connected component 
of Domain(/fc), which takes all values from — oo to oo. We have 



Domain(/fc 



i\ — ^ 

[0,vr]-{^},.ii for odd k 
[0, vr] — {^^^^^} ? for even k 
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Hence the number of real Solutions p of ( |12D with < p < vr is given by 

- 2 for odd k and 



cos^l < ^ 
1 for odd k and |cos^| > ^ , 




The total number is 
N 



N 



2 (A^ - 3) + #{A: = 2j - 1 I 1 < j < ^ 



k-K. N - 

cos — > 

- N 



-}, 



and it is equal to the number of real Solutions {Ai, A2} of BAE for / = 2, including those with the 
value 00, which is the Solutions corresponding to p = ^ for the equation (|12D. Hence the number 
of finite real Bethe Solutions {Ai, A2} for Z = 2 is equal to 



iV, , , , , iV , /cvr, iV-2, , , 

_(Ar_3)-(iV-l) + #{Ä: = 2j-l I l<j<- , |cos— I >^^} . (13) 

Since the total number of complex Solutions of (^) is N , one obtains the contribution of non-real 
Bethe Solutions for / = 2 whose number is given by 

1 for odd k and |cos^| < , 

for odd k and |cos^| > ^ , 

1 for even k and /c 7^ , 
for = , 
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with the total number 
iV- 



N , k-K, iV-2, 
1 - = 2j - 1 I 1 < i < - , I cos — I > -j^} . 



(14) 



Therefore the number of Bethe Solutions for Z = 2 is the sum of ( p^ and (|T^, which is equal to 
— Ci . This coincides with the counting of Bethe 2-states given in . □ 
Example 3. Bethe Solutions for Z = 3 , = 6. The BAE is given by 



A, — n 



m=l J 



Am + i 



Ar 



(j = 1,2,3) , Aj /AfcGC 



Claim: There are exactly 5 Solutions of the above equation, and each of them is invaraint under 
the map A — A. 

First let us consider a Solution invariant under the sign symmetry: A ^ — A. We may assume 



Aj's take the form 



Ai = -As = A 0) , As 



. 



With the variable z in (11), the corresponding equation of zi is 



z'' - 3z^ + 3z - 1 = . 

The above equation has 6 distinct roots, which includes z = 1. The value A corresponding to each 
of the other 5 Solutions gives a Bethe Solution Symmetrie under the change of sign. They contribute 
5 independent states in the Hilbert space Hg, which is a 64-dimensional vector space . On the other 
hand, by Example 1 and 2 of this section together with (iv) of Proposition 1, the total number of 
Bethe states ^'Ar(Ai, . . . , A;) for < / < 2 is equal to 

1x7 + 5x5 + 9x3 = 59. 

Therefore there is no other Bethe 3-state except the Symmetrie ones we have described. Then the 
conclusion follows immeadiately. □ 

Explicit Bethe Solutions of a higher l are difficult to obtain in general for a finite size A^. For 
the rest of this paper, we shall consider only the case 



l = g := — and Xj £ IR for j 
The equation we are going to discuss is the following form: 

Am + i 



1, 



,9 • 



(^^3 + 2 ^Af 



m=l J 



Ar 



, Aj G R {j = l,...,g) 



(15) 



The Bethe vector corresponding to the above Bethe roots leads to the ground state of antiferro- 
magetic Hxxx (i-e. J < ) in — > oo. It is more convenient to consider the logarithmic form of 
the the above equation. By using the relation 



1 X-i 
-logT-— 

i X + 1 



2 arctan(A) + vr , | arctan(A)| < 



TT 



for A e M, 



arctan A 
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we obtain the following equation from (|l5|): 



Q 1 ^ 

arctan 2\j = ^r-^ + ^ ^ arctan(Aj — A^) , 1 < j < g 



k=l 



where Qfs are all integers or all half-integers according to g odd or even, and they are bounded 
by the number Q determined by the relation: 



arctan(oo) = ttIQ + 77) + T7 5Z a'rctan(oo) . 
2 iV 



Hence 



Q 



N_ _ 1 
4 2 



„.._vj — zil -\- J- 1 — 12 n 

and dH) is now equivalent to the equation 

Fj{A) := arctan 2Aj — (vr^j + J2k=i ai"ctan(Aj — A^)) = 



A := 



GM^ , j = l,2,...,g. 



(16) 



Since only a finite number of Solutions can be obtained for BAE (15) by Propostion 1, we obtain 
the following result. 

Lemma 2 . BAE ( |l5|) is equivalent to the equation ( [l6|) , which has at most a finite number of 
Solutions. □ 



Now we are going to show the existence of real Solutions for (16). Define the endomorphism F 
and the linear Involution E of by 



F : 1R9 — >TR3 ^ F(A) : = 



E-.M^ — > IR^ , E{A) 



Denote the £^-invariant vector in M^: 



( 



E{A)j 



V 



, E{A)j :- -Xg+i-j . 



^1 



It is easy to see that F has the following symmetry properties: 
Lemma 3. 

(i) 

F{-A) = -F{A) - 2Tiq for A G R^' . 

(ii) F o E = E o F. □ 

We are going to show the existence of Solutions of the equation (jl^) by the fixed point theory. 
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Proposition 3 . For a sufficient large cube C in 

C = {A = (Ai, . . . , XgY : -a<Xj <a, l<j<g} 

there exists a Solution for the equation F{A) = in A G C. 
Proof. Since 

TT TT 

lim (arctan2x ) = — 

a;-^oo 4 4 

there is a positive number a such that 

TT TT TT 

arctan 2a — 7rq„ = arctan 2a 1 > — . 

^ 4 2iV 4 

Let C be the cube in with 

a > a . 

Denote Cj~ , C~ the faces of C : 

C+ = {A G C I A,- = a} , er = {A G C I Aj = -a} , i = 1, . . . , <? • 
By the inequalities 

TT i ■r-^ , , 9"^ TT 

arctan 2a — irqj > arctan 2q — nqg > — , — arctan(a — Afcj) < — — = — 

A:=l 

one has 

Fj{A) > for A G . 

For A G C~ and j > [^], we have 

qj>0 , -A G C+ 

hence by (i) of Lemma 3, 

Fj{A) = -F,{-A)-27rqj<0 . 

For A G C~ and j < [^], we have 

E{A)€C+^,_^, F,+i_,(ü;(A))>0 

hence by (ii) of Lemma 3, 

F^{A) = -Fg+,_,{E{A))<0. 

Therefore we obtain 

Fj(A)<0 forAGC". 

Thus by Poincare-Miranda fixed point theorem, (see e.g. Q p.p. 12), the conclusion of the propo- 
sition follows immediately. □ 

Remark. (i) One can require the symmetry property on the Solutions in the above proposition. 
Indeed there is a Solution of F{A) = in the intersect of cube C with the hypersurface H, 

H ■.= {A£(S3 \ £;(A) = A} . 

In fact, the map F sends H into itself. Since Aj, 1 < j < [^], form a coordinate System of H, the 
conclusion in the above proof implies that the map 

-^rest ■ H nC — > H 
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also satisfies the the conditions of Poincare-Miranda fixed point theorem, hence it follows the result. 

(ii) The argument given in the proposition also provides an appriori estimate on the location 
of roots of -F(A) = 0, which lies in the following region 

{A I bj <Xj<aj, l<j<g} 

where ^ ^ 

aj = - tan(- + ngj) bj = - tan(— + ngj) . 

□ 

As a corollary of Proposition 3 and Lemma 2 , we have the following result: 
Proposition 4 . There is a Solution of the equation (|l5|). □ 

The uniqueness for the Solution of (|T5| ) should be expected by the thermodynamic nature of the 
Solutions. Let us look a few cases of small N. Note that if {Xj}j^i satisfies the equation (16), they 
satisfy the following equality: 

9 

arctan(2Aj) = . 

For N = A, the above symmetry relation enables one to obtain the Solution of the equation ( [l^ ) 
which is equivalent to 

I Ai + A2 = 

1 arctan(2Ai) = ^ + \ arctan(Ai — A2) 

Hence 



For a general A'^, one can conclude that all the Aj's can not be of the same sign. Indeed one can 
determine signs of Ai and A^ from the first and the last relations in (|l6|): 



-1 1 , ,A^ , TT , . , A l , ,N , TT 

\ ) + ( 1) = , arctan(2A„) > tt( ) - ( 1) 

4 2Ar^ ^ 2 ^2Ar ' y 9J \^ 2N' ^2 ' 2N 



hence 

Ai < < Ag . (17) 
It appears to be the case that certain symmetry properties exist among Aj's, e.g. 

Ai < A2 < • • • < Ap , Aj + Ag+i_j = , 

but the mathematical derivation from the equation (^) seems a difficult problem, even in the case 
of = 6: 

arctan(2Ai) = + g(arctan(Ai — A2) + arctan(Ai — A3)) 
arctan(2A2) = |(arctan(A2 — Ai) + arctan(A2 — A3)) (18) 
arctan(2A3) = | + |(arctan(A3 — Ai) + arctan(A3 — A2)) 

However in the above case, by the analysis of Example 3 in this section, the symmetry property 
for the Solutions holds: 

Ai + A3 = , A2 = . 

Hence Ai satisfies the equation 

5arctan(2A) — arctan(A) = vr . 
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Since the left hand side is a strictly inceasing function of A, there exists an unique Solution of 
Ai, hence the uniqueness of the equation (18). For a larger size A^, the mathematical structure 



of the equation (|16| ) becomes more comphcated that no effective mean could be found at this 
moment for the uniqueness problem. In the next section, we shaU present an plausible, but not 
a mathematically rigorous, argument on the unique ground State Solution for a large but finite N 
based on the thermodynamic limit procedure. 

5 Ground State for Antiferromagnetic XXX Model 



The ground State of the Hamiltonian -f^xXX antiferromagnetic case is the state for the 

JV 

Solution of (|16|). In the thermodynamic limit, one assumes there is a real Solution {Aj} for the 



asymtotic equation of (|16|): 



JV 

arctan2Aj ^ irgj + — arctan(Aj — Xk) , 1 < J < — • 
k=i 

The Qj^s are considered as quantum numbers of the ground state. As N tends to oo, the continuous 
Version of the above relation is obtained by the following Substitution : 

-1 1 

qj — >xe ( — , -) , Xj — > X{x) , 

here A(x) is a monotonic increasing function with A(^) = — oo and A(^) = oo. The density of the 
ground state is now defined to be 

dx 1 
^^^^^ = dÄ ^= nJ1%^xN{X,+,-X,) 

The logarithmic BAE for the ground state now becomes 



arctan2A(x) = nx + arctan(A(x) — X{y))dy. 



Differentiating the above equation with respect to x, one can derive the integral equation for p^: 

vrp.(A) + {-^ * p.)(A) = (19) 
here the convolution of function / and g is defined by 

/•OD 

{f*9)W= / f{X-i,)g{p)di,. 



Then the density, energy , momentum and spin of the ground state in the thermodynamic limit for 
the antiferromagnetic XXX model are given by 

1 Ntt 

P^(^) = ö TTTT' = «^^108 2 , = — (mod 27r), = . 

zcosh(7rA) 2 

( For the details, see || @ ). 

Now we explain the reason for the uniqueness of the ground state for a large but finite N. 
Suppose that {Aj}^^ is a set of Solutions of ( [l^ ) such that the density Pv{X) of the continuous 
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limit A(x) is given by 2cosh(7rA) • Without loss of generality, we may assume there is another set of 
real Solutions {XjYj^i of the equation (16) for a size A^. Then for some fij, one has the expansion: 

-11 N 

^3=^j + J^l^i + 0{j^) ,J = l,...,y. 

We have 

2 1 1-1 

arctan(2A,) + ^^^^^ = vr,, + - g[arctan(A, - A,) + , ^ (,^. _ + 0{-,) , 

hence 

^ . N/2 1 _ 1 

In the continuous limit, 

Xj — > A(x) , — > nix) for X G ^) , 
the above equation becomes 

2fi{x) fl fi{x) - ß{y) 



l + 4A(x)2 l + (A(x) -A(y))2 

Therefore 



dy . 



2/i(A) /i(A)-/xH . w 



1 + 4A2 y_oo 1 + (A - z.)2 

and 



By (19), one has 



M>~)PJ.>~) = -I TT#^'*'- 

, 1 + (A - i/)^ 



By taking Fourier transform, the above equation becomes 



-1 

1+"Ä2' 



which implies p,{x) = 0. Therefore \j and Xj agree up to the order of Repeating the same 
procedure inductively to higher order terms of , one arrives the conclusion that Xj coincides with 
Xj for all j. 
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